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Abstract 

We discuss the heavy-quark contribution to deep inelastic scattering in the scheme with rif = 3,4,5 fixed flavors. 
Based on the recent ABM11 PDF analysis of world data for deep-inelastic scattering and fixed-target data for the 
Drell-Yan process with the running-mass definition for heavy quarks we show that fixed flavor number scheme is 
sufficient for describing the deep-inelastic-scattering data in the entire kinematic range. We compare with other PDF 
sets and comment on the implications for measuring the strong coupling constant a s (Mz). 
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Account of the heavy-quark contribution to deep- 
inelastic scattering (DIS) is an important issue for the 
PDF fits. The charm- and bottom-quark contributions, 
which are mostly relevant for the existing data analy- 
sis, are substantial at small values of the Bjorken vari- 
able x forming an important contributions to the HERA 
data on the inclusive structure functions (SFs). Further- 
more, the semi-inclusive SFs, which correspond to the 
subprocess with the heavy quarks in the final state, can 
provide an additional constraint on the small-x gluon 
distribution [1 1. However, the calculation of the heavy- 
quark production cross section is hampered within per- 
turbative QCD because the higher-order corrections are 
extremely involved for the case of two scales, the lep- 
ton momentum transfer Q 2 and the heavy-quark mass. 
At present the QCD corrections to the massive Wilson 
coefficients are known up to the NLO only J2] |3] |H. 
The partial NNLO corrections stemming from the soft- 
gluon threshold re-summation have been also calcu- 
lated for the neutral-current (NC) @ and charged- 
current (CC) heavy -quark production (6). For the case 
of electro-production they are numerically important at 



small x and Q 2 reaching O(10%) in the kinematic re- 
gion of HERA. The NNLO threshold re-summation cal- 
culations of Ref. [5 1 were recently updated [7, 8| and 
now they include all threshold-enhanced logs and the 
Coulomb term. The heavy-quark mass appearing in 
the massive Wilson coefficients of Refs. d [2] [4] [8] 
corresponds to the pole mass which emerges in the 
QCD Lagrangian. The pole-mass definition provides 
a straightforward way for the perturbative calculations, 
however, it is not ideal for phenomenology since the 
pole mass is quite sensitive to the QCD radiative correc- 
tions. This shortcoming is eliminated when the heavy- 
quark mass is defined in the MS -scheme, similarly to 
the strong coupling constant a s . The heavy quark MS - 
masses are conventionally parametrized at the scales of 
the heavy-quark mass itself. These scales are close to 
the typical DIS hard-scattering scale therefore their per- 
turbative stability is greatly enhanced if compared to 
the pole-mass definition. The massive Wilson coeffi- 
cients re-calculated in terms of the running mass also 
demonstrate improved perturbative stability and the re- 
duced renormalization/factorization scale sensitivity at 
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Figure 1: The CC inclusive cross section data obtained at HERA 1 10 1 with the positron beam (left) and the electron beam (right) 
compared with the running-mass 3-flavor scheme predictions based on the NNLO PDFs of Ref. |9| (solid lines). The same 
predictions without the c-quark contribution taken into account are given for comparison (dashes). 



the HERA kinematics (9). In this range only the charm- 
and bottom-quark production are relevant. The c- and 
/7-quark MS -masses averaged over different determina- 
tions are 

m c (m c ) = 1.27 ±0.08 GeV (1) 

and 

m b (m b ) = 4.19 ±0.13 GeV, (2) 

respectively ifTTl . The uncertainties in these values are 
essentially smaller than the spread in the heavy-quark 
masses values employed in the various global PDF fits, 
cf. Tab. 1 in |9|. Therefore the running -mass definition 
provides a footing for the consolidation of those PDFs. 

Taking advantage of the running-mass definition 
and new massive NNLO corrections we update our 
ABKM09 PDFs, which were obtained within the fixed- 
flavor-number (FFN) scheme at NLO and NNLO IfTHl . 
and produce in this way a new PDF set, ABM1 1. Apart 
from this theoretical improvement we add to the fit 
recently published HERA data. Firstly, the separate 
HI 02) and ZEUS OH data sets on the NC inclu- 
sive SFs are replaced by much more accurate combined 
HERA data ifTUl . The unprecedentedly small uncer- 
tainty of 1-2%, including the normalization error of 
0.5%, achieved for the combined HERA inclusive data 
allows to reach better constraints on the small-* PDFs 



as compared to the ABKM09 analysis. Only the NC 
combined HERA data with Q 2 < 1000 GeV 2 are used 
in order to exclude contributions due to the Z-boson 
exchange without loosing the statistical significance of 
the NC HERA data for our fit. Besides, we include 
into the fit the CC HERA data obtained for the elec- 
tron and positron beams |10|. The CC HERA sample 
ranges up to Q 2 = 15000 GeV 2 . To demonstrate that 
the treatment of the charm production within the FFN 
scheme is applicable at that large transverse momentum 
we compared the predictions based on the variant of the 
NNLO ABKM09 fit performed with the running-mass 
scheme 0. The value of m c (m c ) = 1.18 ± 0.06 GeV 
is obtained in this fit with the constraint of Eq. ([T]i im- 
posed. Taking this value of m c (m c ) we find very good 
agreement with the CC HERA data in the whole range 
of Q 2 , cf. Fig. [T] At small x the charm contribution to 
the inclusive CC cross sections is quite significant, sim- 
ilarly to the case of the NC DIS. Since the CC charm 
production is initiated by the strange quarks mainly, this 
sample provides an additional constraint on the small-x 
strange sea distribution, which is at the moment defined 
by the neutrino-nucleon DIS data only |fl5l . In the same 
way the charged-lepton initiated CC data can be used to 
determine the value of the c-quark mass, in particular 
employing the potential of the planned high-luminosity 
EIC facility [16|. Finally, we add to the fit the inclusive 
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HI data obtained in a special HERA run with high in- 
elasticity y achieved [17]. At the kinematics probed in 
this run the inclusive SFs are sensitive to the structure 
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Figure 2: The pulls versus x for the HI NC inclusive DIS c.s. 
data of Ref. 1 17 1 binned on the momentum transfer Q 2 in units 
of GeV 2 w.r.t. the ABM1 1 NNLO fit. The data points with dif- 
ferent inelasticity y still may overlap in the plot. The inner bars 
show statistical errors in data and the outer bars the statistical 
and systematical errors combined in quadrature 1181 . 

function Fl at small x and Q 2 and thereby shed light on 
the details of the small-x QCD dynamics, in particular 
the shape of the gluon distribution at small x. 

Similar to the ABKM09 analysis, in ABM 11 the 
HERA data are supplemented by the ones on dimuon 
production in the (anti)neutrino-nucleon DIS and the 
Drell-Yan process in order to separate the quark PDFs 
by flavors. Also we employ the fixed-target inclusive 
DIS data obtained in the NMC ||24), the BCDMS j25] 
1261 . and the SLAC experiments l27ll28ll29ll30ll3Tll32l. 
The fixed-target data allow to constrain the PDFs at 
large x. However, modeling of the higher-twist terms 
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Figure 3: The predictions for the semi-inclusive structure func- 
tion F L 2 C at different values of the momentum transfer Q 2 ver- 
sus x for the ABM 11 PDFs in the NNLO (solid curves), the 
ABM1 1 PDFs in the NLO (dashes), and the NN21 PDFs (H) 
in the NLO (dashed dots), all taken in the 3-flavor running- 
mass scheme with the the value of m c = 1.27 GeV II II . The 
NNLO predictions based on the MSTW (20) and JR09 ED 
PDFs (not displayed in the plot) strongly overlap with the 
ABM1 1 ones. The HI data of Ref. 1221 given in the plot are 
extracted with HVQDIS code of Ref. |23) . 



is required since they contribute to the DIS SFs at small 
Q 2 651 . The corrections for nuclear effects fl34l are 
also taken into account in the analysis of fixed target 
data when relevant, cf. [18] for details. The overall 
quality of the fit is quite good. For 3036 data points 
used in the fit the value of^ 2 is 3391 and 3378 for the 
NLO and NNLO variants, respectively. For the com- 
bined HERA data [10) including the NC and CC sam- 
ples the value of x 2 /NDP = 537/486 is obtained at 
NNLO with full account of the systematics error corre- 
lations. The pulls of NC HERA data w.r.t. the fit do not 
demonstrate a statistically significant trend versus Q 2 . 
This gives additional justification for application of the 
FFN scheme to the analysis of the DIS data at realistic 
kinematics, cf. also [35 1. At small x and Q 2 the per- 
turbative QCD corrections rise and they are particularly 
big for the NNLO corrections to the massless Wilson 
coefficients for Fl [36). With account of these correc- 
tions the high-y HI data 11171 . which are quite sensitive 
to the contribution from Fl are described pretty well, 
with the value of x 2 = 137 for 130 data points. Also the 
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HI data do not suggest a violation of the conventional 
QCD evolution down to x ~ 1CT 5 , cf. Fig. [5] In this 
way we do not confirm a hint on the evolution kernel 
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Figure 4: The data on F L versus x obtained by the HI col- 
laboration 1 17 1 confronted with the 3-flavor scheme NNLO 
predictions based on the different PDFs (solid line: this anal- 
ysis, dashes: JR09 ED, dots: MSTW08 (20)). The NLO 
predictions based on the 3-flavor NN21 PDFs 1191 are given 
for comparison (dashed dots). The value of Q 2 for the data 
points and the curves in the plot rises with x in the range of 
1.5 -=-45 GeV 2 GU- 

re-summation effects observed in the NNPDF fit ll38l . It 
is worth noting in this connection that the interpretation 
of the low-x DIS inclusive data is sensitive to the treat- 
ment of the heavy-quark contribution. While we em- 
ploy in the analysis of DIS data the FFN scheme with 3 
flavors, in many other PDF fits l20ll37l[39l the variable- 
flavor-number (VFN) scheme is used. The VFN ap- 
proach to the description of the DIS is based on the 
asymptotic expressions for massive Wilson coefficients 
taken in the 3-flavor FFN scheme Il40l I4H . which are 
valid at Q 2 » w?, where m/, is the heavy-quark mass. 
In this limit the power corrections vanish and the struc- 
ture functions come down to convolutions of the mass- 
less Wilson coefficients with the PDFs, which now are 
defined in the 4-flavor scheme including also the heavy- 
quark PDFs. The matching conditions between 3- and 
4-flavor schemes are known up to the second order in 
a s Pol |4D and the first third-order corrections are also 
available ll42l [43). However, in order to employ the 
VFN scheme in the analysis of realistic data, which 
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Figure 5: The 3-flavor gluon distributions obtained in the vari- 
ous NNLO fits at the scale of fi = 1.2 GeV (solid line: this 
analysis, dashes: JR09 ED, dots: MSTW08 EO), dashed 
dots: NN21 t37l ). 



range down to Q 2 ~ m 2 h , one has to take into account 
the power corrections as well. This task is conceptually 
difficult within the VFN formalism therefore the power 
corrections are commonly modeled using the combina- 
tion of the massive FFN and massless VFN Wilson coef- 
ficients Il44) . The models which are obtained in this way 
pretend to describe the full-range kinematics of existing 
DIS data; therefore they should reproduce the 3-flavor 
scheme results at small Q 2 . A particular shape of the 
massive VFN coefficient functions at Q 2 ^f> m? is sub- 
ject to a particular choice and there are numerous VFN 
scheme prescriptions, which differ basically by the de- 
gree of smoothness provided for matching with the FFN 
scheme. In cases the smoothness is achieved by intro- 
ducing empirical parameter(s), which control dumping 
of the massless VFN term at small Q 2 l45ll . It is worth 
noting that the uncertainty in a particular choice of these 
parameters is propagated into the uncertainties in PDFs 
obtained with such prescriptions. Moreover, additional 
parameters, which appear in the coefficient functions of 
Ref. 1 45 1 do not enter into the QCD anomalous dimen- 
sions. Therefore the QCD factorization may be broken 
in this way. In contrary, the NNLO S-ACOT-^- prescrip- 
tion ll46ll is explicitly based on the factorization theo- 
rem and does not contain damping factors. However the 
matching smoothness is not guaranteed in this case and 
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the S-ACOT-^- based calculations overshoot the FFN 
scheme results at Q 2 ~ ml, cf. Fig. 4 in |46| . On the 
other hand, the BMSN prescription |40| ensures smooth 
matching with the FFN scheme at Q 2 = m 2 h inherently, 
without additional parameters. Furthermore the results 
of the PDF-analysis are not sensitive to the choice of ei- 
ther the VFN scheme as the BMSN prescription or the 
FFN scheme fl2l . 

To a certain extent different VFN scheme prescrip- 
tions can be sorted out with the data on the semi- 
inclusive structure function F?£, which correspond to 
the DIS sub-process with two charmed quarks in the fi- 
nal state. Strictly speaking this observable is infrared 
unsafe due to the non-singlet term contributing to per- 
turbative QCD calculations starting from the NLO. Rig- 
orous infrared safety of F| c can be restored only with 
an additional soft cut, e.g. imposed on the cc invariant 
mass |47|. At the same time the non-singlet contribu- 
tion in F"' is numerically small and it can be safely dis- 
regarded in the analysis of the existing DIS data. The 
recent HI data on F£ c Il22l are compared to the pre- 
dictions based on different PDFs in Fig. [3] To provide 
a comparison consistently we employ the FFN scheme 
with the running-mass definition, taking the correspond- 
ing 3-flavor PDFs provided by different groups and the 
world-average value of m c (m c ) = 1.27 GeV [QQ |. This 
value was also used in our analysis, however, the data 
of Ref. [22 1 are not included into our fit. The calcula- 
tions are performed with the open-source code OPEN- 
QCDRAD Q which contains the implementation of the 
DIS Wilson coefficients up to the NNLO. The NNLO 
ABM 11 predictions are in a good agreement with the 
data. Also they agree with the combined HI and ZEUS 



data on F| c , which range wider in kinematics |48 1. The 
NLO ABM 11 calculations are very close to the NNLO 
ones that reflect the improved perturbative stability of 
the running-mass definition scheme. The NNLO pre- 
dictions based on the MSTW flU and JR09 (211 PDFs 
are in a good agreement with ours and with the HI data, 
while the NLO NN21 lfl9l predictions systematically 
undershoot the data. The NN2 1 PDFs are obtained with 
the FONLL prescription of the VFN scheme |49|. The 
FONLL variant of the VFN scheme is conceptually sim- 
ilar to the prescription of Ref. [45 1, however, it allows 
less flexible modeling of the massive coefficient func- 
tions at Q 2 ^> m 2 h . In any case the VFN modeling should 
reproduce the FFN scheme results at small Q 2 . How- 
ever, the trend demonstrated by the NN21 predictions is 
opposite and they also diverge from the HI data at small 
Q 2 - 

'http : //www-zeuthen . desy . de/~alekhin/OPENQCDRAD 



The high-y HI data ifTTll also discriminate different 
PDFs. The predictions for Fi at low x computed with 
various NLO and NNLO PDFs are compared with the 
HI data in Fig. |4] As for the comparison with F| c 
above, we employ the FFN scheme with running-mass 
definition and the publicly available 3-flavor PDFs. The 
NNLO ABM 11 and JR09 predictions obtained in this 
way demonstrate good agreement with the data in the 
whole range of x, while the MSTW predictions under- 
shoot the data at low x. This is correlated with the strong 
fall-off of the MSTW gluon distributions at small x tak- 
ing negative values at x < 10~ 3 at small scales, cf. 
Fig. [5] The observed discrepancy cannot be attributed 
to the impact of a particular choice of the VFN scheme 
prescription in the MSTW fit since the low-x tail of the 
HI data correspond to small Q 2 , where the VFN scheme 
employed in the MSTW fit reproduces the FFN one. 
Due to this the small-x MSTW gluon distribution should 
move up and consolidate with the JR09 and ABM ones 
once the HI data flTTj are included into the MSTW fit. 
The same is also valid for the NN21 PDFs, which over- 
shoot the HI data at small x and undershoot them at 
larger x, cf. Fig. [5] 

The value of a s is a necessary ingredient of all QCD 
calculations and it should be consistent with the PDFs 
employed. In our analysis we provide this consistency 
fitting a s simultaneously with the other parameters of 
the data model. In this way we obtain the value of 



a s (M z ) = 0.1134 ±0.0011 (at NNLO) (3) 



and 



a s (M z ) = 0.1180 ±0.0012 (at NLO) 



(4) 



with the uncertainties corresponding to the 68%C.L. 
They are calculated with the standard statistical crite- 
rion Ax 2 = 1 taking into account correlations of the 
systematic uncertainties in the data within the covari- 
ance matrix approach ll53l . Both values in Eqs. ( 3]4 l 
are at variance due to the (9(0.005) scale variation error 
at NLO. Once a s is fitted simultaneously with the PDFs 
and the higher-twist terms in the DIS SFs the errors in 
Eqs. ( 3|4 i take into account uncertainties in the latter as 
well. The ABM1 1 value of a s is in a very good agree- 
ment with our earlier ABKM09 value and the error is 
somewhat smaller due to the improved accuracy of the 
HERA data. The profiles ofy 2 versus a s , obtained in the 
NLO and NNLO variants of the ABM11 analysis with 
the value of a s fixed and all other parameters fitted, are 
nearly parabolic at the minimum with the shapes deter- 
mined by Eqs. ( |3|4| >, cf. Fig. [8] This gives an additional 
justifications of using the standard statistical criteria to 
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Figure 6: The charm- (left) and the bottom-quark (right) PDFs obtained in the global fit: The dotted (red) lines denote the ±lcr band 
of relative uncertainties (in percent) and the solid (red) line indicates the central prediction resulting from the fit with the running 
masses |9|. For comparison the shaded (gray) area represents the results of ABKM09 1 12] 1181 . 
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Figure 7: The data on charged-lepton asymmetry versus the lepton pseudorapidity 77 obtained by the ATLAS |50| (left panel) and 
CMS |51 1 (right panel) experiments compared to the NNLO predictions based on the DYNNLO code of Ref. j52| and the ABM1 1 
NNLO PDFs (shaded area showing the integration uncertainties). The ABKM09 NNLO predictions are given for comparison by 
dashes, without the integration uncertainties shown [18]. 



estimate the parameter uncertainties in our fit. The x 2 are also nearly parabolic, cf. Fig. [9] The value of a s 
profiles for the separate data sets employed in the fit preferred by the HERA and BCDMS data are in a good 
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agreement, while the NMC and SLAC data prefer some- 
what smaller and bigger value, respectively. The SLAC 
and the NMC data are sensitive to the higher-twist con- 
tribution, cf. discussion in ||T8"ll54"l . Note, the cut on the 
hadronic invariant mass W, which is commonly used in 
the global PDF fits, does not allow to get rid of the im- 
pact of the higher-twist terms on a s . Indeed, in a variant 
of our fit with the cut of W 2 > 12.5 GeV 2 imposed and 
all higher- twist terms fixed at we obtain a s (M z ) = 
0.1191 + 0.0006 at NNLO. This is substantially bigger 
than the value of a s in Eq. Q. It is worth noting that in 
this way we approach the values of a s obtained in the 
PDF fits [20 55], performed with no higher-twist terms 
are taken into account. The value of a s is also sensi- 
tive to the treatment of the correlated uncertainties in the 
data. E.g. in the NNLO MSTW fit [20 1 the HERA and 
NMC data prefer value of a s (M z ) > 0.12, contrary to 
our findings, cf. Fig. [9] Note that in [20 1 the NMC and 
HERA systematics errors are combineid in quadrature. 
To study impact of this approximation we performed a 
trial NNLO ABKM09 fit with the same treatment of the 
NMC and HERA systematics and found that the value 
of a s (Mz) shifted up by +0.0029 as compared to the 
nominal ABKM value. Many other aspects of our anal- 
ysis, which may affect the value of a s , are also different 
from [20, 55): basic relations for the DIS cross sections, 
data normalization, etc, cf. Ref. [ 18 1 for a detailed dis- 
cussion. These differences make a detailed comparison 
of our results with (20] [55 ] difficult. 

While the 3-flavor FFN scheme is nicely sufficient 
for description of the DIS data, the 5-flavor factoriza- 
tion schemes is most often justified for the collider phe- 
nomenology in view of much bigger scales involved. In 
cases the 4-flavor scheme may be also relevant. The 
4(5)-flavor ABM 11 PDFs are matched with the 3(4)- 
flavor ones at the scale of m c and m/,, respectively, em- 
ploying massive OMEs taken in the running-mass def- 
inition (9|- The 4- and 5-flavor PDFs at larger scales 
are generated by means of the massless QCD evolution 
with these boundary conditions. The heavy-quark PDFs 
are particularly sensitive to the values of m c , mi,. Taking 
advantage of the MS definition we fix them at the PDG 
values of Eqs.(T]2 1. Due to the change in the heavy fla- 
vor treatment and impact of the new data included into 
the fit, the ABM 11 heavy-quark PDFs differ from the 
ABKM09 ones, cf. Fig. [6] The uncertainties in the 
heavy-quark PDFs are to a large extend defined by the 
uncertainties in m c , Mb- In our fit the latter are calcu- 
lated as a sensitivity of the fitted data supplemented by 
the PDG constraints of Eqs.(T]2i to m c ,mt,. The uncer- 
tainty in nib obtained in this way coincides with one in 
Eq.([T]i, while for m c it reduces to 0.06 GeV due to im- 
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pact of the inclusive HERA data. The uncertainties in 
the heavy-quark PDFs estimated with these constraints 
on m c ,m/, are essentially reduced as compared to the 
ABKM09 case. This improvement is especially impor- 
tant for the hadronic single-top production driven by 
the initial-state b-quarks and for the Higgs production 
through the vector-boson-fusion channel, which is sen- 
sitive to the c-quark distribution |56|. 

The hadronic jet production provides an additional 
constraint on the PDFs, in particular on the large-x 
gluon distribution ll20l [39l \l9l . However, the calcula- 
tion of the full NNLO QCD corrections to this process 
is still in progress (see If57]l58l and references therein). 
This precludes a consistent use of the Tevatron jet data 
in our NNLO PDF fit. Nevertheless, in order to check 
any potential impact of the jet Tevatron data on our 
PDFs we have performed trial variants of the NNLO 
ABKM09 fit with the Tevatron jet data added (59). The 
NLO QCD corrections EH |6T) and the partial (soft 
gluon enhanced) NNLO corrections due to threshold re- 
summation 1 62] have been computed with the FastNLO 
tool 116311641 . 

In general, the Tevatron jet data overshoot the 
ABKM09 predictions, nevertheless they can be 
smoothly accommodated in the fit. The typical value of 
X 2 INDP sb 1 is achieved with account of the error corre- 
lations for the jet data sets of [ 67 68 65 . 66 1 once they 
are included into the NNLO ABKM09 fit. Meanwhile 
the various data sets demonstrate a somewhat different 
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Figure 9: The ^-profile versus the value of a s (Mz), for 
the HERA data CH [TTJ, the NMC data (24), the BCDMS 
data HSlUa, and the SLAC data t27ll28ll29l[30ll3T1[32l. all 
obtained in variants of the present analysis with the value of 
a s fixed and all other parameters fitted (solid lines: NNLO fit, 
dashes: NLO one) (TH). 



Experiment 


a s (M z ) 




NLO 


NNLO* 


DO 1 jet 


0.1190 ±0.0011 


0.1149(12) 


DO 2 jet 


0.1174(9) 


0.1145(9) 


CDF 1 jet (cone) 


0.1181(9)9 


0.1134(9) 


CDF 1 jet (JfcJ 


0.1181(10) 


0.1143(9) 


ABM 11 


0.1180(12) 


0.1134(11) 



Table 1: The values of a s (M z ) based on including individual 
data sets of Tevatron jet data |65 66 67 68 1 into the analy- 
sis at NLO. The NNLO* fit refers to the NNLO analysis of 
the DIS and DY data together with the NLO and soft gluon re- 
summation corrections (next-to-leading logarithmic accuracy) 
for the 1 jet inclusive data, cf. 



trend with respect to the ABKM09 predictions. E.g., the 
off-set of the DO inclusive jet data [67] does not depend 
on the jet energy Ej and therefore may be attributed to 
the impact of the currently missing full NNLO correc- 



tions, cf. Figs. 10 and 11 In contrast, for the CDF data 
of [65 1 obtained with the k T jet algorithm the pulls rise 
with Ej and can be reduced only by means of a modifi- 
cation of the PDF shapes. 

The values of a s extracted from the trial ABKM09 
fits with the Tevatron jet data included are compared 



with the nominal ABKM09 value in Tab. [T] At most, 
they are bigger by l<x, while for the CDF cone jet algo- 
rithm data |66| the central value of a s is even the same. 
The predictions for the light Higgs production cross 
section, which are defined by the gluon distribution at 
x < 0.1, are also not very sensitive to the constraints 
coming from the Tevatron data, cf. Tab. [2] The impact 
of the Tevatron jet data on the large-x gluon distribution 
is more significant. However, in this context we note 
that the trend of the first LHC data on the jet produc- 
tion with respect to the various PDF predictions is dif- 
ferent from the Tevatron measurements. The ABKM09 
predictions are in better agreement with the CMS and 
ATLAS inclusive data of l70l ITTI than the predictions 
based on the PDFs of E0ll39l[T9l , which were tuned to 
the Tevatron inclusive jet data. Jet data from LHC is still 
subject to large systematic errors, though. Note also that 
the Tevatron dijet and 3-jet production data are in good 
agreement with the ABKM09 predictions ll72l . in con- 
trast to the case of inclusive jet production at Tevatron. 
These ambiguities in the data as well as the limitations 
in the current theoretical treatment prevent the use of 
hadronic jet data in our fit. 

The data on W/Z productions being produced by the 
LHC experiments also may help to improve the PDF 
accuracy. The charged-lepton asymmetry data ll50l I5D 
obtained by the ATLAS and CMS experiments are com- 
pared to the NNLO predictions based on the ABM11 
PDFs in Fig. [7] All differential distributions for W- and 
Z-boson production are computed with the fully exclu- 
sive NNLO program DYNNL0 [52 1, which allows to take 
into account the kinematical cuts imposed in the exper- 
iments (cf. Fig. [7J. The overall agreement with both 
experiments is sufficiently good. At values of rj ~ 1.5 
for the lepton pseudo-rapidity rj the data show a differ- 
ent trend with respect to the predictions, however the 
discrepancy is within the data uncertainties. Prelimi- 
nary data on the charge-lepton asymmetry at large ra- 
pidities obtained by the LHCb collaboration ll73l are 
also in good agreement with the ABKM09 predictions. 
To check the impact of the LHC charged-lepton asym- 
metry data on our fit, we have performed a variant of 
the ABM 11 analysis which consists of adding the data 
of l50l ED to the fit. We have found that the impact of 
those data is only marginal in view of still big experi- 
mental uncertainties. 

Summary. We have discussed different schemes for 
the treatment of the heavy flavor component in DIS. 
In the kinematic region Q 2 ^ all schemes need to 
match to the 3-flavor scheme, which allows for a con- 
sistent comparison of the various PDF sets. We have 
performed this benchmark exercise with the help to the 
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cr{H)[pb] 


ABKM09 


DO 1-jet inc. 


DO di-jet 


CDF 1-jet inc. 
(cone) 


CDF 1-jet inc. 

(h) 


Tevatron(1.96) 
LHC(7) 


0.770(50) 
14.34(41) 


0.859(29) 
14.68(29) 


0.833(27) 
14.69(27) 


0.815(25) 
14.11(28) 


0.842(25) 
14.44(27) 



Table 2: The predicted cross sections for Higgs boson production in ggF with m H = 120 GeV at Tevatron ( = 1-96 TeV) and at 
LHC ( -\fs = 7 TeV) from NNLO variants of the ABKM09 fit 1121 corresponding to Tab. [T] The uncertainty in brackets refers to 
the lcr standard deviation for the combined uncertainty on the PDFs and the value of a s (M z ). The values in bold correspond to the 
published result |69|. 



OPENQCDRAD code for standardized precision com- 
parison. Within the framework of the ABM11 PDF 
analysis, which uses the running mass scheme for the 
heavy quarks, we find that the FFN is completely suffi- 
cient for describing the existing DIS data. 

We have also given detailed information on theoret- 
ical and experimental improvements which are of rele- 
vance for the low-x PDFs, especially for the gluon and 
we have briefly sketched the implications for LHC phe- 
nomenology, so that differences of ABM1 1 with respect 
to other PDF sets can be explained. More benchmark 
comparisons with the help of OPENQCDRAD can be 
found elsewhere 1 18 1 and are expected in the future. 
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